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ABSTRACT 

We present high-precision timing of five millisecond pulsars (MSPs) carried out for more than 
seven years; four pulsars are in binary systems and one is isolated. We are able to measure the 
pulsars' proper motions and derive an estimate for their space velocities. The measured two- 
dimensional velocities are in the range 70-210 km s"^, consistent with those measured for other 
MSPs. We also use all the available proper motion information for isolated and binary MSPs to 
update the known velocity distribution for these populations. As found by earlier works, we find 
that the velocity distribution of binary and isolated MSPs are indistinguishable with the current 
data. Four of the pulsars in our observing program are highly recycled with low-mass white 
dwarf companions and we are able to derive accurate binary parameters for these systems. For 
three of these binary systems we are able to place initial constraints on the pulsar masses with 
best-fit values in the range 1.0-1.6 Mq. The implications of the results presented here to our 
understanding of binary pulsar evolution are discussed. The updated parameters for the binary 
systems studied here, together with recently discovered similar systems, allowed us to update 
previous limits on the the violation of the strong equivalence principle through the parameter |A| 
to 4.6x10"^ (95% confidence) and the violation of Lorentz-invariance/momentum-conservation 
through the parameter Idsj to 5.5xl0~^° (95% confidence). 

Subject headings: binaries: close — stars: neutron — stars: pulsar — relativity 



1. Introduction 

Pulsars are believed to be born with spin pe- 
riods of ~0.1 s and gradually slow down as they 
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age due to the loss of rotational kinetic energy 
in the form of electromagnetic radiation. These 
"normal" pulsars make up the bulk of the ob- 
served popula tion and ^1,900 of the m are cur- 
rently known ( Manchester et al. 2005f ^. On the 
other hand, a group of old, fast-spinning pul- 
sars is observed ('^200 pulsars). It is believed 
that these old pulsars are formed from the trans- 
fer of mass and angular momentum from a pre- 
vious or present c ompanion in a binary system 
(|Alpar et al.lll982l ). These pulsars are generally 
seen as "recycled" members of the population: old 
pulsars that have been spun-up and brought back 
to an active, pulse-emitting life thanks to their in- 
teraction within a binary system. The fastest spin- 



^See also the ATNF Pulsar Catalogue: 
|http: //www. atnf . csiro .au/research/pulsar/psrcat/l 
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ning pulsars known have spin periods of Pg < 0.01 
s, so-called "millisecond pulsars" (MSPs), and are 
thought to have been produced in this manner. 
The measured characteristics of the members of 
these binary systems and their orbital parameters 
provide valuable insights into the formation and 
evolution of these systems. See, e.g., Stairs (2004) 
and Lorimer (2008) for general reviews of binary 
pulsars and their scientific importance. 

The measurement of a pulsar's proper motion 
can be used to estimate its space veloci ty (e.g., 
Hobbs et al.1l2005t IChatteriee et al.ll2009l ). Such 
measurements are important for a variety of sci- 
entific questions, including estimating the distri- 
bution of natal kicks imparted to proto-neutron 
stars by the supernova (SN) explosion that cre- 
ated them. A variety of mechanisms have been 
propo sed that can give rise to these natal kicks 
(e.g.. Spruit fc PhinnevI 1998t Kusenko fc Segre 
199a Ijessner et al.l |2005I )~ Another key ques- 
tion is whether isolated MSPs have a similar ve- 
locity distr ibution to those still in binary sys- 



Roche-lobe overflo w phase (|van den Heuvellll994 : 
Tauris et al.1l2000> ). 



tems (e.g.. iTauris fc Bailes 11996'; Toscano et alJ 
19991: iLommen et al.l Ipood : iHobbs et al. 200^ 
McLaughlin et al.l l2005l )~ In general, we expect 
MSPs to have lower system velocities than the rest 
of the pulsar population since, after the SN explo- 
sion, the binary system must have remained intact 
to spin up the neutron star. However, for isolated 
MSPs, the companion must eventually leave the 
system or be evaporated. Arguments for both 
lower and higher velocities for isolated MSPs have 
been given in the literature ( e.g., To scano et al.l 
19991 : iMcLaughlin et al]l2005[ ). Given that only 
a small number of isolated MSPs with measured 
proper motions are known (<10 objects), addi- 
tional measurements are very important. 

Many subclasses o f pulsar binary systems ar e 
now recognized (e.g.. Stairs 2004t Lorimer 20081 ). 
The broader distinction made is between those 
pulsars with high- mass companions (e.g., another 
neutron star) and those with lower mass compan- 
ions (e.g., white dwarfs - WDs). In the case of pul- 
sars with WD companions, various subgroups are 
generally identified. For example, mildly-recycled 
pulsars (Ps ^ tens of milliseconds) in a tight or- 
bit (orbital periods < a few days) with high- 
mass WD companions (m2 ^ 1 Mq) are thought 
to arise from a common-envelope evolution or 
from periods of ultra-high mass transfer during a 



A more straightforward evolution is thought to 
apply for MSPs in long orbits {Pi, > 4 days) with 
low-mass WD companions (m2 < 0.3 Mq), gen- 
erally call ed wide-orbit binary millisecond pulsars 
(WB MSP; [Rappaport et al.lll995l : lTauris fc Savoniie 
19991 ). Here, as the companion evolves and over- 



flows its Roche-lobe during the red-giant phase, 
mass spirals onto the neutron star and forms an 
accretion disk. A stable, long-lived mass trans- 
fer phase is expected to take place, producing 
nearly circular o r bits ( eccentricities of e ^ 10^^- 
10"^). IPhinnevI (|l992l ) predicted that these sys- 
tems should exhibit an orbital period-eccentricity 
{Pb — e) relationship based on the expectation 
that convective eddies in the envelope of the red 
giant will produce nonzero values of the eccen- 
tricity. The point at which mass transfer ceases 
and e freezes depends on the size of the red giant 
envelope, which will in turn determine the size of 
the orbit and thus the orbital period. In addition, 
the size of the envelope is also thought to be re- 
lated to the mass of the red giant's core, which 
eventually contracts to form a WD. Therefore, an 
orbital period-core mass (Pf, — m2) relationship is 



also expected in these systein s (jRappaport et al 
199,4 iTauris fc Savoniielfl999l) . 



The WBMSP systems also provide impor- 
tant tests for theories of gravity. For exam- 
ple, the strong equivalence principle (SEP) states 
that all neutral test masses fall with the same 
acceleration in an external gravitation field, 
i.e., it states that the gravitational and iner- 
tial masses of self-gravitating bodies are identi- 
cal {mg/rrii = 1). Binary pulsar systems allow 
us to test fo r SEP violations in the hmit of high 



self-gravity C Damour fc Schafed 119911 : IWexl 11997 



2000l : IStairs et &L.2Q(^ if the binary components 



experience different accelerations in the gravita- 
tional field of the Galaxy, a forced eccentricity is 
imparted to the system along the projected direc- 
tion of the external force onto the orbital plane. 
Binary pulsars with low companion mass, small 
eccentricity and long orbital periods are ideal for 
SEP violation tests. 

Another important test of gravitational theories 
involves the post-Newtonian parameter as which 
is associated with the violation of momentum con- 
servation and the existence of preferred frames 
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(Lorentz invariance; I Willi 119931 ) . In general rel- 
ativity (GR), as = 0. The most observable ef- 
fect of a possible deviation from this GR predic- 
tion is thought be a non-zero self-acceleration for 
a rotating body in a direction perpendicular to its 
spin axis and perpendicular to its velocity with re- 
spect to the absolute rest frame (jBell fc Damour 
19961 ). In the case of binary systems, each compo- 
nent will experience self-acceleration. These self- 
accelerations perturb the orbital dynamics, lead- 
ing to a forced eccentricity and polarization of the 
orbit along a fixed direction. 

Here we present results obtained from long- 
term timing of 5 MSPs, four of which are WBM- 
SPs (PSR J1853-t-1303, PSR J1910+1256, PSR 
B1953-t-29, PSR J2016-f 1948) and one of which 
is isolated (PSR J1905-t-0400). We use our re- 
sults to study the space velocities of MSPs, binary 
evolution models and equivalence principle tests. 
In 31] we describe the observations performed and 
data analysis carried out. In ^we present our im- 
proved timing solution for these pulsars, including 
the measured proper motions. In §3.11 we discuss 
the implications for the velocity distribution of iso- 
lated and binary MSPs and discuss their implica- 
tions for binary evolution models of MSPs. In H'6.2\ 
we derive constraints on the component masses for 
three systems and we discuss these results in light 
of evolution models. In we present updated 
upper limits of equivalence principle violations us- 
ing WBMSPs. Finally, in fgl we summarize our 
findings and point to future directions of this re- 
search. 

2. Observations and Data Analysis 

We have conducted high-precision timing on 
four WBMSPs and one isolated MSP. We collected 
data from two observatories and a total of five data 
acquisition systems. Here we describe the observ- 
ing setups used at each telescope. A summary of 
the observations for each pulsar is given in Table 

m 

2.1. Arecibo 

All pulsars were observed with the 305-m 
Arecibo telescope in Puerto Rico. Th e Wideband 
Arecibo Pulsar Processors (WAPPs; iDowd et al 



cept PSR J2016-hl948 for which all four WAPPs 
were used in some observations. They were oper- 
ated in online folding mode with 32 /is sampling 
and 192 lags near 1400 MHz and 96 lags at 2700 
MHz (for PSR J2016+1948 a sampling time of 
128 /zs and 128 lags were used) . The Arecibo Sig- 
nal Processor (ASP: [Pemorest 2007t ) was used for 
all pulsars except PSR J2016-hl948. ASP provides 
0.25/iS complex sampling in two orthogonal polar- 
izations. These data were coherently de-dispersed 
in software using 16 or 24 frequency channels, each 
with 4 MHz bandwidth. The polarisations were 
later summed and the signal folded at the pulsar 
period. The ASP observations were flux calibrated 
with a pulsed noise diode of known strength and, 
when available, with observations of a standard 
flux calibration source. A typical observing ses- 
sion involved collecting multiple integrations of 1 
or 3 min long on each pulsar with ASP and the 
WAPPs for a total of up to 30 min of data. The 
data from the short integrations were aligned and 
summed using a preliminary timing model. For 
the WAPPs, all data were summed to obtain at 
single profile for each observation. For ASP, a 
separate profile was obtained at each frequency 
channel for each observation. 

For PSR J2Q16-hl948 we also used d ata from 
the P enn State Pulsar Machine (PSPM; ICadwelll 
19971 ). an analogue filterbank with 128 x 60 kHz 



200fl[ ) were used to observe all the pulsars. Three 



frequency channels. The power level for each 128 
channels was sampled every 80 fis and stored to 
tape. The data were subsequently folded and 
aligned multiple times as the ephcmcris for the 
pulsar was being refined. The PSPM and WAPP 
data for this pulsar were summed every few min- 
utes and multiple profiles were obtained for most 
epochs. 

To provide a long time baseline and improve 
the measured proper motion for PSR B1953-f29 
we also used d ata taken with the Mark II system 
()Rawlevl[l986 '). a dual-polarization 32 x 30 kHz 
filterbank spectrometer. Here, the outputs from 
opposite polarizations were summed and sampled 
at the pulsar period. Data were averaged over in- 
tervals of 1-2 minutes and stored for off-line pro- 
cessing. A total of around 1 hour of data were 
collected at each epoch and the data summed to 
obtain one profile for each observatiord. 



of the four WAPPs were used for most pulsars, ex- 



2 Additional data for PSR B1953+29 from Arecibo (MJD 
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2.2. Parkes 

We used the 64-m Parkes telescope in NSW, 
Australia to observe pulsars PSR J1853+1303, 
PSR J1905+0400 and PSR J1910+1256. These 
observations were taken using the Parkes ana- 
logue filterbank centered at 1390 MHz with 512 
X 0.5-MHz frequency channels sampled every 0.25 
ms. Two polarizations were recorded and summed 
in hardware for each frequency channel. The data 
were subsequently folded off-line using a prelim- 
inary ephemeris and summed to obtain a single 
profile for each observation. 



A time of arrival (TOA) was found for each 
observation by cross-correlating the profil es with 
a hig h signal-to-noise standard template ([Taylor 
19921 ). The recorded observatory times were cor- 
rected to UTC time by using data from GPS satel - 
lites. The JPL DE405 ephemeris (|Standishll20"ol 
was used for barycentric corrections. The soft- 
ware package TEMPCd was used to find the tim- 
ing solution for each pulsar by including astromet- 
ric, binary and spin parameters as needed to ar- 
rive at a phase-connected solution (where every 
rotation of the star over the span of the obser- 
vations is accounted for). In order to fit for any 
instrumental or standard template profile differ- 
ences, we fit for arbitrary time offsets between 
each instrument (for ASP we have also allowed 
for time offsets between each 4 MHz channel to 
fit for any profile changes across its bandwidth). 
A change in dispersion measure (DM) over time 
was detected for PSR J2016-I-1948 and marginally 
for PSR J1905-I-0400 and PSR B1953-I-29 (see Ta- 
ble mi . Finally, the measured uncertainties were 
scaled by a small telescope-dependent amount to 
ensure a timing fit with xt 



3. Results and Discussion 

In Figure [T] we show the standard pulse pro- 
files for each pulsar at 1400 MHz. Table [D shows 
the timing solutions derived from our work and 
Figure [2] shows the timing residuals derived from 
these solutions. The timing solutions successfully 



49129—49255) and the Effelsberg telescope in Germany 
(MJD 49768-50460) were also available but did not add 
significantly to the results and were left out of our analy- 
sis. 

^http:/ /tempo. sourceforge.net/. 



model the measured TO As and leave no significant 
trends in the residuals. Pulsars PSR J1853-I-1303 
and PSR J1910-t-1256 have very low root-niean- 
square (rms) values and are now part of a long- 
term timing program to detect and study gravita- 
tional waves using an array of well-timed puls ars 
( Hobbs et al.ll2010t iDemorest et al ] l2009ll201l[ ). 

PSR J1853-f 1303 and PSR J1910-t-1256 were 
discovere d by the Parkes multib eam pulsar sur- 
vey (e.g., Manchester et al. 200ll) an d their tim- 
i ng so lutions were first reported by IStairs et al. 
([20051). The longer data span possible with the 



Parkes data and the high quality of the Arecibo 
data allowed us to improve the timing solutions for 
these pulsars (especially for their binary parame- 
ters) and, for the first time, report a measurement 
of their proper motions. For PSR J1910-I-1256, 
we are also able to measure a secular change in 
its projected semimajor axis (a similar, though 
marginal, measurement was also made for PSR 
J1853+1303). 

PSR B1953-f29 was discovered while per- 
forming a systematic search for radio pulsars 
using Arecibo in position error boxes from 7- 
ray sources reported by the COS B satellite 
(|Boriakoff et all 11983). Pr evious timing solu- 



tions for PSR B195 3-I-29 have been repo r ted b y 
Rawlev et aP (|l988l ) and IWolszczan et al.l (|2000l ). 



Here we have been able to use Arecibo data span- 
ning 25 yrs to derive a much improved timing 
solution for this pulsar, and especially so for its 
proper motion measurement. PSR J1905-I-0400 
was also dis covered by the Par kes multibeam pul- 
sar survey (jHobbs et al.l 120041 ) and is one of only 
<20 known isolated MSPs. Here we are able to 
measure for the first time a proper motion for this 
pulsar. 

PSR J2016-I-1948 was discovered in t he Arecibo 
430 M Hz Intermediate Latitude Survey ([Navarro et al 



20031 ). The discovery dataset for this pulsar cov- 



ered about a year (taken in 1999) and were enough 
to determine that the pulsar was in a binary 
system with an unusually long orbital period of 
635 days. However, deriving a complete timing 
solution from this initial dataset was not possi- 
ble. It was later discovered that the original 1999 
data have large systematics, most likely the re- 
sult of being folded with an inacurate estimate 
for the spin period. These data have been left 
out of our analysis and all subsequent data were 
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taken in observing modes that allow for re-folding 
and re-aligning as the pulsar ephemeris was being 
improved. The current timing solution leaves no 
systematics in the derived TOAs and has been 
correctly predicting new TOAs for many years. 
We are therefore confident that we have found the 
most precise timing solution currently available for 
this system. Our timing solution in Table [2] shows 
that the pulsar is indeed recycled and in a nearly 
circular orbit, likely the result of mass transfer and 
tides as its companion was going through its giant 
phase. We are also able to measure the proper 
motion for this system and a secular change in the 
projected semimajor axis. 

None of the pulsars show a significant value 
of annual parallax (see Table [2]). Only PSR 
J1853-I-1303 has a marginal parallax measurement 
with a large error. It is possible that further ob- 
servations of this pulsar with improved timing pre- 
cision will be able to produce a more constraining 
parallax measurement. In addition, none of the 
binary pulsars show a measurable Shapiro delay. 
The residuals obtained from the best-fit timing so- 
lution are shown as a function of binary phase in 
Figure |31 

3.1. Millisecond Pulsar Velocities 

The high precision obtained by our timing 
study allowed us to measure a statistically sig- 
nificant value for the proper motions of all five 
pulsars (see Table [2]). We have used these mea- 
surements to study the velocity distribution of 
MSPs, both isolated and those in binary sys- 
tems. The pulsar population in general has been 
found to have large space velocities with a mean of 
^300 km s-i (iLvne k Lorimeij|l994l : iHobbs et af 



20051 ). Recycled MSPs appear to be on the low 
end of the velocity distribution, with a mean of 
~90 km s~^. In addition, no significant difference 
has been found between the velocities distribu- 
tions of isolated MSPs and those still in binary 
systems (jHobbs et al.l l2005t iLommen et al.l 120061 ) 



despite the fact that an additional evolutionary 
stage (the disruption of the binary) has occurred 
in the former. 

Now we revisit the velocity distribution of 
MSPs, which we take to be those with periods 
of P < 0.01 sec and are therefore fully recycled. 
From the new timing solutions presented in Table 
El only PSR J2016+1948 is not a fully recycled 



MSP and its implied 2D velocity of 96 km s ^ (at 
the implied DM distance of 2.5 kpc) was not in- 
cluded in the following analysis. PSR J1905-I-0400 
studied here is particularly important, as it is one 
of only ten isolated MSPs with measured proper 
motions. Our sample then consists of 10 isolated 
MSPs and 27 binary MSPs for a total of 37 pul- 
sars. We have combined the measured proper 
motions with the available distance estimates to 
calculate the pulsars' 2D space velocities in their 
respective local standard of rest at the pulsar lo- 
cation. To do this we have corrected for Solar 
motion and used a peculiar velocity for the Sun 
of Fq=13.4 km s-i |Dehnen fc Binnevll 19981 ). We 



have also assumed a flat Galactic rotation curve 
with a galactocentric distance for the Sun of Rq=8 
kpc and a Galactic rotation velocity of 222 km s~ ^ 
(|Eisenhauer et aLlliool iDehnen fc Binnevlll998[ ). 
A flat Galactic rotation curve is thought to be a 
good approximation for distances from the Galac- 
tic centre of >3 kpc (Oiling & Merrifield 1988; 
the pulsars we used have distances to the Galactic 
center of >4 kpc). The resulting 2D velocities 
in the pulsars' standard of rest after correcting 
for Solar and Galactic motion, ¥20, are shown in 
Tablets 

Most pulsars have distance estimates from tim- 
ing measurements of their dispersion measure 
(DM) combined with a model of the free elec- 
tron distribution in the Galaxy ( Cordes fc Lazij 
I2OO2I) . In general, distances derived using this 
method are thought to have a ~25% error, im- 
plying a minimum similar error on the estimated 
velocities. For individual pulsars, the distance er- 
ror could be much larger than that. The errors on 
the estimated pulsar velocities shown in Table |3] 
were derived using Monte Carlo simulations with 
10,000 runs per pulsar. For these simulations, 
pulsar distances were drawn from Gaussian dis- 
tributions centered on the values listed in Table 
|3] with a width of 25% of the central valu£0 (for 
pulsars where non-DM distances are available the 
corresponding distance errors were used). Pul- 
sar proper motions were then drawn using Gaus- 
sian distributions with central values and widths 
derived using the values listed in Table [H In 
practice, the largest error contribution to the esti- 



^This uncertainty is consistent with the ~20% estimate of 
uncertainty in distances due to unmodeled inhomogeneities 



m 



the interstellar medium model l|Cordes fc Lazi* 
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mated pulsar velocities are the associated distance 
errors. We then used the velocities in Table [3] to 
study the velocity distribution of MSPs. 

Figure S] shows the normalized histograms of 
the 2D velocities of binary (solid line) and isolated 
MSPs (bold dotted line) in our sample. The av- 
erage velocities are found to be 108±15, 113±17 
and 93±20 km s'^ for all MSPs, binary MSPs 
and isolated MSPs, respectivel}lf|. The average 
2D velocities without correcting to the pulsars' 
standard of rest are 88±12, 96±15 and 68±16 
km s~^ for all, binary, and isolated MSPs, re- 
spectively. The upd ated velocities are consistent 
with previous work: iHobbs et all (|2005[ ) reported 
2D uncorrected velocities for binary and isolated 
MSP s of 89ibl5 and 76±16 k m s~\ respectively 



with Lommen et al.l (|2006l ) and iMcLaughlin et al 



(|2005l ) reporting similar values. We then find that 
the average velocities of binary and isolated MSPs 
are consistent with each other. To test whether 
the two samples are consistent with arising from 
the same distribut ion we use t he Kolmogorov- 
Smirnov (KS) test (|Massevlll95lP . A KS test of 
the two corrected velocities results in a probability 
of 62% that they are drawn from the same distri- 
bution. For the uncorrected 2D velocity measure- 
ments, the two distributions have a KS probability 
of 75% that they are drawn from the same distri- 
bution. 

We therefore conclude that there is no statis- 
tically significant difference between the velocity 
distributions of isolated and binary MSPs with the 
current statistics. However, we also note that due 
to selection effects our sample is biased towards 
nearby, low-velocity pulsars. It is therefore possi- 
ble that the lack of difference between the velocity 
distribution of isolated and binary MSPs is due to 
our observing the low velocity tail of these distri- 
butions, which in reality could be quite different. 
Higher number statistics (particularly for isolated 
MSPs) will allow for a more detailed study of such 
effects in the future. In addition, more distant 
MSPs are now being discovered in current surveys 



^The errors in the average velocities shown in this section 
represent the standard errors of the mean. 

^ While the KS test doesn't take the error estimates into 
account, it is one of the most useful and general methods for 
comparing two samples. Detailed simulations to account 
for the errors are beyond the scope of this paper and can 
be carried out in future work. 



(e.g., in addition to PSR J1903-I-0327 in Table [3 
the velocities for two distant MSPs will also be 
published by Deneva et al. 2011). While measur- 
ing the proper motions of distant objects will most 
likely require large amounts of telescope time, they 
represent significant additions to our sample. 
Furthermore, we note that Tauris fc Bailed 



(|1996[ ) presented the expected space velocities of 
binary MSPs using various evolutionary models. 
In their simulations, binaries with shorter peri- 
ods tend to have larger velocities since the final 
velocity of the system depends on the separation 
of the components at the time of the supernova 
explosion. However, this correlation is fairly weak 
and asymmetries in t he explosion would easily 



wipe out this effect. iHobbs et al] (|2005l) found 



no evidence for a correlation between the velocity 
of binary MSPs and their binary period^ We 
now briefly revisit this idea and in Figure [5] we 
plot the binary periods versus implied velocity for 
the binary pulsars listed in Table [3] This figure 
should be compared t o the plots in Figure 2 of 
Tauris fc BailesI (|l996[ ). For system with Pb < 2 



days and Ph > 2 days we find average 2D velocities 
of 135±52 km s~^ and 107±14 km s^^, respec- 
tively (uncorrected 2D velocities have averages of 
120±45 km s^^ and 85±11 km s'^). 

We therefore find no significant difference in the 
velocities of short- and large-period binary MSPs. 
However, we caution that only a handful of the 
former systems are known. While the very large 
velocity of PSR B1957-I-20 could be explained by 
these models using asymetries in the supernova 
explosion, they would have a particularly hard 
time explaining the small implied velocities for 
PSR J0751+1807 and PSR 2051-0827 given their 
very short orbital per iods. In addition, we did 
simulations using the iTauris fc BailesI ()l996l ) bi- 
nary period— velocity relationship taking into ac- 
count the effect of random projections towards us 
of these velocities. We find that the large scat- 
ter in the relationship and the random projections 
of these velocities would most likely wash out any 
effect. Therefore, while we find no evidence that 
the relationship is present (and this might indeed 



iHobbs et all ||2005| ') used a definition of P < 0.1 s and 
P < lO"^'^ s s~^ for recycled pulsars, thus including a mix 
of companion types in their sample. Here we use P < 0.01 s, 
resulting in binaries with mostly helium white dwarf com- 
panions. 
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be very difficult to achieve, even if it exists), at 
the same time we cannot rule it out. It is clear 
that additional work is needed to understand the 
evolution of binary MSPs. Obtaining additional 
velocity measurements for these pulsars will help 
to constrain evolutionary models. 

3.2. Component masses and change in 
projected semimajor axis 

For PSR J1910+1256 and PSR J2016+1948 we 
were able to measure a change in the projected 
semimajor axis, x=dx/dt (see Table [2]). For PSR 
J1853+1303, the measurement of x was marginal 
and will be discussed at the end of this section. For 
PSR B1953+29 only an upper limit was measured. 
Here we define x=aismi/c, where ai is the semi- 
major axis of the pulsar orbit, i is the inclination 
angle of the angular momentum vector of the or- 
bit relative to the Earth-pulsar line of sight (LOS), 
and c is the speed of light. The measured values 
for x are -1.8(5) x lO"!'^ and 8.3(14) x 10"" for 
PSR J1910+1256 and PSR J2016-hl948, respec- 
tively. In principle, a non-zero value for x could 
arise from a change in ai , i or a combination of the 
two. However, we argue that the measured values 
likely arise due to the pulsars' high proper motion 
inducing a change in our LOS to these binaries. 

In the case that a change in ai is being ob- 
served, GR predicts a value of |a'i| ^ 5x10^^"^ 
and -10-24 for PSR J1910- j-1256 and PSR 
J2016+1948, respectively (see IPetersI Il964 for 
the required expression). These values are many 
orders of magnitude below the observed value 
of i;. In addition, for typical binary astrophys- 
ical processes, a non-zero value for |a'i/a| is ex- 
pected to have a similar order of magnitude as 
\Pb/Pb\- No significant value for i\ was found 
for any of our pulsars, but allowing for a mea- 
surement in our timing solution results in a value 
of P6=-2(4)xl0-" and -1(2) xlO'^ s s^^ for 
PSR J1910-M256 and PSR J2016-H1948, respec- 
tively. Using the values for Pb listed in Table [2] 
we find limits of \Pb/Pb\ <1.2xl0"" s^^ and 
<5xl0-" s-i for PSR J1910-hl256 and PSR 
J2016-I-1948, respectively. Again, these values are 
a few orders of magnitude lower than expected 
from the measured values of x. 

We therefore propose that the observed val- 
ues of X must arise from apparent changes in the 
orbital parameters due to the proper motion of 



the b inaries ( Arzoumanian et al. 19961 Kopeikin 



19961 ). The strength of this geometric effect can 



be calculated using: 



1.54x10 ^^x^ cot i sin 6* ss 



(1) 



where x is the projected semimajor axis in units 
of seconds, is the total proper motion of the sys- 
tem in units of mas/yr and 9 is the unknown angle 
between the position angle of the proper motion 
and the position angle of the ascending node of the 
pulsar's orbit. The measured values of x can then 
be used to constrain the inclina t ion an gle of the 
system i. Following iNice et al.l (|200ll ). we have 



used a Monte Carlo simulation to constrain the 
values of i that satisfy Equation [T] for both PSR 
J1910+1256 and PSR J2016-f 1948, resulting in la 
values of 44°(36°-52°) and 36°(27°-45°), respec- 
tively. 

In addition, numerical studies of the evolution 
of neutron star binaries in long-period orbits with 
low-mass white dwarf companions point to a rela- 
tionship between the final orbit al period, Pb, and 



the mass of the companion, (jRappaport et al 



199,4 iTauris fc Savoniiel[l999h . An overall agree- 



ment with these results has been found in the 
available data, although these relationships ap- 
pear to overestimate m2 for systems with long pe- 
riods and provide conflicting results for systems 
with short periods (jStairs et al.l 120051 ). 

Keeping these caveats in n iind, we used the 
Ph — W 2 relationship found bv iTauris &: Savoniid 
(|l999f ). together with the inclination angle con- 
strains from i, to provide constraints on the 
masses of the PSR J1910+1256 and PSR J2016+1948 
systems. The derived constraints on the com- 
panion masses are TO2 = 0.30-0.34 M0 and 
0.43-0.47 Mq for PSR J1910-I-1256 and PSR 
J2016-I-1948, respectively. These values and those 
implied by the mass functions are shown in Figure 
m Restricting 7712 to lie in the values implied by 
the Pb — m2 relationship and using the inclination 
angles from x we find la values for the pulsar 
masses of wi = 1.6±0.6 Mq and 1.0±0.5 Mq for 
PSR J1910+1256 and PSR J2016-I-1948, respec- 
tively. 

For PSR J1853+1303, the measured value of 
X is marginal (see Table [2]) but can still be used 
to derive initial constraints on the system masses. 
Following the same procedure as outlined above 
results in a 1(7 confidence interval for the incli- 
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nation angle of this system of 48°(33°-58°). In 
addition, the Ph — m2 relationship produces com- 
panion masses of m2 = 0.33-0.37 M©. Combining 
these results, we derive la values for the mass of 
PSR J1853+1303 of mi = 1.4±0.7 Mg. 

The derived mi values are not very constrain- 
ing, though fully within the expected mass ranges 
for neutron stars. Given that PSR J2016-t-1948 is 
only one of three WBMSP with A > 200 days (see 
Table H]) , further constraining the masses of this 
system by independent measurements and contin- 
ued timing can provide a valuable constrain to bi- 
nary pulsar evolution models. 

3.3. Theories of Gravity: Tests 

We have used the improved timing parameters 
for the four WBMSPs studied here, in addition 
to recently discovered systems, to update impor- 
tant tests of GR and other theories of gravity. In 
particular, we have modeled the forced eccentric- 
ity that would be imparted on the binary systems 
due to violations of the strong equivalence princi- 
ple (SEP) using the parameter A, and the forced 
eccentricity that would be imparted due to viola- 
tions of Lorentz invariance/momentum conserva- 
tion using the parameter a^. In GR, both param- 
eters are predicted to be identically zero and ds is 
also predicted to be zero in most other theories of 
gravity 

For the SEP test parameter A, the additional, 
forced eccentricity imparted on the binary orbit 
is exp ected to be of the form (jPamour fc Schafer 
19911) : 



by (|Bell fc Damouij|l996l) 



eF,A 



A- 



(2) 



2FGM{2T:/Pb)^ 

where c is the speed of light, F is unity in GR 
and a function of mi and m2 in alternate theories, 
G is Newton's constant in GR, M=mi+m2, Pb 
is the binary period and |g^| is the projection of 
the Galactic acceleration vector onto the orbital 
plane at the location of the pulsar. Here, the total 
observed eccentricity is then predicted to be eo^s = 
e^r-l-eir A, with the "natural" eccentricity e^r and 
the angle 6 between e^r and e^^^^A being additional, 
unknown parameters. 

For the Lorentz invariance/momentum conser- 
vation test parameter da, the forced eccentricity 
added to the binary orbit is expected to be given 



"3- 



247r P GM 



sin/3 



(3) 



where cp~~2E^°'^ /mic^ is the gravitation self- 
energy fraction of the pulsar (the so-called "com- 
pactness", taken to be approx i mately 0.21mi: 
[Pamour fc Esposito-Fares3ll992t iBell fc Damour 
1 19961 ). B is the unknown angle between the pul- 
sar system's absolute velocity V (with respect 
to the reference frame of the cosmic microwave 
background) and the pulsar's spin vector. 

We have used the above expressions and a 
Bayesian analysis to derive probability distri- 
butions for A and da given the measured bi- 
nary/pulsar parameters and additional estimates 
for the remaining unknown param eters. The 
proced ure was described in detail in IStairs et al.l 
( 2005 ) and we summarize it her^. For each pulsar 
j, we find the probability density functions (pdf) 
p(|A||iDj,/) andp(d3|-Dj, /) for probable values of 
A and da given each pulsar's data Dj (see Table 
m) and prior relevant information /. 

For example, for the SEP parameter A we can 
write for each pulsar: 

pi\A\,dj\D,,I)<^p{D,\\Ald,,I)xp{\Ald,\I) 

(4) 

where dj represents the relevant parameters for 
this test, namely i, m2, fl, d, cn, and 9. For 
these parameters we perform Monte Carlo simula- 
tions when they are not directly measured or con- 
strained through timing or other methods. For m2 
we use twice the range given by the P h — m2 re- 
lationship of Tauris fc Savoniig ( 19991) . For cosi 



we assume a uniform distribution between and 
1, and combine this value with the measure mass 
function and m2 to provide a value for the pulsar 
mass mi (only systems with mi values between 
1.0 and 2.5 are keptjl. For Q we use a uniform 
distribution between 0° and 360°. For d we use 
a Gaussian cen t ered o n the DM estimate using 
Cordes fc Lazio (|2002h and assuming an average 



uncertainty of 25%, or a Gaussian centered on the 
parallax measurement if available. The integrals 



*We have also fixed some small bugs in the IStairs et al.l 
i2005 ) code that did not significantly alTect their results. 

^For the pulsars presented here with new measured values 
for X, we have not included the implied orbital constaints 
as they have large error bars. 
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over the remaining unknown parameters e n and 9 
are computed separately using the measured val- 
ues of Bobs and u! and the implied values of eii-^A- 
A pdf for the parameter was similarly derived; 
for this we also need estimates for the 3D velocity 
of each pulsar and used Gaussian distributions in 
each dimension centered on the Galactic rotational 
velocity vector at the pulsar loca tion with widths 



of 80 km s (jLvne et al.l 119981 ) or, when avail- 



able, we use Gaussian distributions for the proper 
motions to get the transverse velocities. For |A|, 
the parameter space 10~^ < |A| <0.1 was sam- 
pled uniformly in steps of 2x10"^ and for jdal the 



-19 . 



parameter space 10 < jasj <5xl0 was sam- 
pled uniformly in steps of 1x10^^^. 

Binary systems suitable for these studies are 
required to have large periods and small eccen- 
tricities so that additional relativistic effects are 
negligible. Large values of /e and P^/Pe have 
therefore been used as a general selection c har- 
acter i stic for choosing appropriate systems (jWexl 
2OOOI : Stairs et al. 2005 ). In addition, the systems 
must be old enough and have large enough u so 
that the orientation of their orbits can be assumed 
to be random and that the projection of the Galac- 
tic acceleration vector on the orbit can be assumed 
to ha ve been constant over the lifetime of the sys- 
tems ( Damour fc Schafeij[l99ll : IWexlll997[ ). While 
some pulsars might individually provide low limits 
for these tests, we use all 27 available systems in 
order to provide a more conservative upper limit 
that incorporates the assumptions made on the 
population as a whole. 

Currently, a total of 27 WBMSPs are available 
to test the above effects and their properties are 
listed in Table HI The pdfs for each pulsar are 
shown in Figure [7] for the SEP parameter |A| and 
in Figure [8] for the Lorentz invariance/momcntum 
conservation parameter Idaj. Since each pulsar 
represents an independent test of these parame- 
ters, we can multiply the individual pdfs to obtain 
a total pdf from our sample of pulsars. 

Using all the systems in Table [H for |A| we 
derive a 95% upper limit of 4.6x10""^, which 
represents a 20% impr ovement from the value 
derived by IStairs et al.l ((2005). Two new pul- 
sars are particularly constraining for this test: 
PSR J1711-4322 and PSR J1933-6211, which 
together with the improved parameters for PSR 
J1853-t-1303 have significantly contributed to the 



reduced upper limit for |A| (the secondary peak 
at low values of |A| of ~10~'^-^ in the product pdf 
shown in Figure [7] is mainly due to these pulsars). 
For PSR J1711-4322 alone, the 95% upper hmit 
for |A| is 5.6x10^"'. Since pulsars test gravita- 
tional theories in the regime of strong fields, future 
improvements to the above limit are important. 
The fact that two pulsars discovered in the last 
five years were able to significantly contribute to 
this test is encouraging and raises the possibility 
that additional discoveries, and improved param- 
eters for the objects already known (particularly 
PSR J1711-4322 and PSR J1933-6211), will im- 
prove the limit further. 

For Ids I, using the updated sample of pulsars 
we derive a 95% upper limit of 5.5x10^-^°. This 
limit is hi gher than the value of 4xl0~^° de- 
rived bv iStairs et al.l (j2005l ). We believe that the 
higher value better reflects the limits of this tech- 
nique to constrain jdal when a larger sample of 
pulsars is available. The most constraining pul- 
sars for this test currently are PSR J1713-I-0747 
and PSR J1853-I-1303 with very similar limits of 
2.8x10-20 and 3.1x10-^0 (95% confidence), re- 
spectively. The limits derived here are about 13 
orders of magnitude smaller than Solar System 
values (Will 1993) and again test the strong field 
limit. Further discoveries and ongoing study of 
present systems (particularly PSR J1853-I-1303) 
will help to place additional constraints on this 
test of gravitational theories. 

4. Conclusions 

We have presented updated timing solutions for 
five MSPs, four of which are in binary systems 
and one which is isolated. The high precision and 
large time span of the observations used allowed 
us to measure the proper motion of these pulsars. 
The implied 2D space velocities in each pulsar's 
standard of rest lie in the range 70-210 km s~^. 
We studied the available 2D velocities of binary 
and isolated MSPs and find that their velocity dis- 
tributions are indistinguishable with the current 
data. For PSR J1910-hl256 and PSR J2016+1948, 
we are able to measure a significant rate of change 
of the semimajor axis which we attribute to a ge- 
ometrical change in our line of sight to the pulsars 
due to their high space velocities. We are then 
able to put initial constraints on the mass of these 
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pulsars of 1.6±0.6 Mq and 1.0±0.5 M© for PSR 
J1910+1256 and PSR J2016+1948, respectively. 
For PSR J1853+1303 we measured a marginal rate 
of change of the scmimajor axis, resulting in an es- 
timate for the pulsar mass of 1.4±0.7 Mq. 

We are also able to place updated constraints 
on violations of the SEP and Lorentz invari- 
ance/momentum conservation using an updated 
list of binary pulsars in wide orbits with small ec- 
centricities. Using a total of 27 pulsars we derive 
an upper limit for the SEP violation parameter 
|A| of 4.6x10"^ (95% confidence) and an upper 
limits for the Lorentz invariance/momentum con- 
servation violation parameter \as\ of 5.5xl0~^° 
(95% confidence). Further discoveries and ongo- 
ing study of present systems will help to improve 
these limits. 
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J1653+1303 





Fig. 1. — Standard profiles for each pulsar at 1400 MHz obtained by combining all the Arecibo data used 
in our analysis. For PSR J2016+1948 we used the WAPP data and for all other pulsars we used the ASP 
data. The x-axis shows one pulse period and the y-axis shows arbitrary flux units. 
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Fig. 2. — Post-fit timing residuals for each pulsar. From top to bottom we show: PSR J1853-I-1303, PSR 
J1905+0400, PSR J1910+1256, PSR B1953+29 and PSR J2016+1948. For all pulsars, the black TOAs are 
those obtained from ASP (except for PSR J2016+1948, where black data show the TOAs obtained with 
PSPM). In all cases, the red TOAs are from the WAPPs. For all pulsars, the blue TOAs are those obtained 
from Parkes (except for PSR B1953-I-29, where the blue TOAs are those from the Mark 11 instrument). 
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Fig. 3. — Same as for Figure [5] but showing only the binary pulsars in our sample with their residuals plotted 
as a function of orbital phase. 
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Fig. 4. — Normalized histograms of the 2D velocity distribution of binary (solid line) and isolated (dashed 
line) MSPs. The errors for each bin are given by the propagated measurement errors. 
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Fig. 5. — Orbital period versus 2D velocities for binary MSPs. 
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Fig. 6. — Constraints on the inclination angle, i, and companion mass, TO2, for PSR J1910+1256 (top) 
and PSR J2016+1948 {bottom). Solid line are constraints derived from the measured mass function of the 
systems. Vertical dashed lines represent inclination angle constraints derived from the measured x values 
(centre line is the median likelihood and outer lines represent the la likelihood limits). Horiz ontal dot- 
dashed lines are the m2 values derived from the Pb — m2 relationship in Tauris fc Savonii^ ( 19991 ). See 
for details. 
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Table 1: Summary of Observations 



Telescope Instrument 


Number of 


MJD 


Center 


Eflcctive 




TOAs 


Range 


Frequencies 


Bandwidth 








(MHz) 


(MHz) 


PSR J1853+1303: 










Arecibo ASP 


494 


53370-55105 


1400 


64 or 96 




23 


54999-55105 


2350 


64 or 96 


WAFF 


41/32/38 


53061-55134 


1170/1370/1470 


50/50/50 




7/8/6 


54882-55134 


2650/2750/2850 


100/100/100 


Parkcs Filtcrbank 


46 


52606-54023 


1390 


256 


PSR J1905+0400: 










Arecibo ASF 


371 


53370-54808 


1400 


64 or 96 


WAFF 


29/23/39 


52279-54808 


1170/1370/1470 


50/50/50 


Parkcs Filtcrbank 


87 


51492-53835 


1390 


256 


PSR J1910+1256: 










Arecibo ASF 


430 


53370-55105 


1400 


64 or 96 




49 


54882-55105 


2350 


64 or 96 


WAPP 


32/26/22 


53187-55171 


1170/1370/1470 


50/50/50 




8/9/9 


54882-55171 


2650/2750/2850 


100/100/100 


Parkes Filterbank 


48 


52602-54062 


1390 


256 


PSR B1953+29: 










Arecibo ASP 


205 


53912-55105 


1400 


64 or 96 




15 


54967-55106 


2350 


64 or 96 


WAFF 


17/13/13 


53912-55134 


1170/1370/1470 


50/50/50 




2/2/2 


54882-55171 


2650/2750/2850 


100/100/100 


Mark II 


47 


46112-49096 


430 


0.96 


PSR J2016+1948: 










Arecibo PSPM 


324 


52456-53591 


433 


7.68 


WAFF 


775/355/557/199 


52939-55392 


1170/1310/1410/1510 


100/100/100/100 
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Table 2: Measured and Derived Parameters for the Observed Pulsars 

Parameter PSR J1853+1303 PSR J1905+0400 PSR J1910+1256 PSR B1953+29 

Right Ascension (RA), a (J2000.0) 18:53:57.319174(8) 19:05:28.273436(16) 19:10:09.701479(8) 19:55:27.87600(3) 

Declination (Dec), 5 (J2000.0) 13:03:44.0784(2) 04:00:10.8830(6) 12:56:25.5074(3) 29:08:43.4659(5) 

Epocli (MJD) 54000 53700 54000 54500 

Data span (MJD) 52606.1-55134.8 51492.2-54808.8 52602.2-55171.8 46112.6-55134.9 

Proper motion in RA, /iQ=acos(5 (mas yr"^) -1.68(7) -3.80(18) 0.21(10) -0.9(1) 

Proper motion in Dec, Mi = i (mas yr"^) ... -2.94(12) -7.3(4) -7.25(12) -4.1(1) 

Annual parallax, tt (mas) 1.0(6) <2.5" <0.7" <7" 

Spin frequency, !/ (s^^) 244.3913778653740(15) 264.242346143483(16) 200.658805375034(1) 163.04791306911(2) 

Spin frequency derivative, !> (s~^) -5.2060(5) X 10" '^'^ -3.425(1) X 10"^** -3.900(2) X 10" "^"^ -7.901(3) X 10"^** 

Dispersion measure, DM (pc cm"^) 30.5701(6) 25.6923(12) 38.0701(8) 104.501(3) 

First DM derivative, DM (pc cm"^ yr"^) ... <|4xl0"'*|" -l.l(7)xlO"^ < 16x10""!" -4.7(2.5) X lO"^"" 

Orbital period, (days) 115.65378643(2) 58.466742029(12) 117.34909728(4) 

Project semimajor axis, £C (It-s) 40.7695198(3) 21.1291036(3) 31.4126903(8) 

Rate of change of a;, i (Is-s s"^) 1.7(7)xl0"" -1.8(5)xl0"" <4xl0"'-'''' 

Eccentricity, e 2.3691(12) X 10"^ 2.3018(3) X 10"* 3.3025(5) X 10"" 

Longitude of periastron, ui (deg) 346.60(4) 106.014(9) 29.485(8) 

Epoch of periastron. To (MJD) 54046.78(1) 54079.318(1) 54444.267(2) 

Number of TOAs 695 549 633 316 

Weighted rms residual (/isec) 1.5 5.4 1.8 3.8 

Derived Parameters 

Spin period, P (ms) 4.09179744490025(2) 3.78440484691321(7) 4.9835840178364(1) 6.1331666053350(1) 

Spin period derivative, P (s s"i) 8.7163(8) X 10"^^ 4.905(2) X 10"^^ 9.687(4) X 10"^^ 2.9734(1) x 10"^" 

Surface magnetic field, B (G) 1.9x10** 1.4x10** 2.2x10** 4.3x10** 

Spindown luminosity, E (10^^ ergs s"^) 5.1 3.5 3.1 5.1 

DM-dcrived distance, doM (kpc) 2.1 1.7 2.3 4.64 

Characteristic ago, Tc (Gyr) 7.4 12.2 8.1 3.3 

Mass function, /i (Mq) 0.00543963576(7) 0.002962840(12) 0.00241678837(2) 

Minimum companion mass^, m2 (Mq) 0.24 0.19 0.18 

Total proper motion, /J (mas y"^) 3.39(11) 8.24(36) 6.98(14) 4.2(1) 

Galactic angle of proper motion*^, 6^ 274° 270° 210° 197° 

Note— Values in parentheses are uncertainties in the last digits shown, which are twice the formal errors quoted by TEMPO after scaling the TOA 
uncertainties to obtain xi^ —1. Right ascension values are in hours, minutes, and seconds and declination in degrees, arcminutes, and arcseconds. 
For all pulsars, the DE405 ephemeris was used and the recorded observatory ti mes were c orrected to TT(BIPM). 
" Value shown represents a Ax^ ~ 6.6 from best fit, representing a ~3cr limit <IAvnilll976H . 

^ For PSR J2016+1948, a second DM derivative {DM) was also measured with a value of 3.8(1.4) X 10"" pc cm"^ yr"^. A less significant value for 
DM was also measured for PSR B1953+29 giving -2.1(1.6) X 10"" pc cm"^ yr"^. 
^ Assuming a pulsar mass of mi — 1.35Mq. 
Clockwise from Galactic North. 



PSR ,12016+1948 



20:16:57.44349(6) 
19:47:51.5882(12) 

53000.0 
52456.2-55392.2 
1.28(26) 
2.83(34) 
<4.5'' 
15.3987376281305(13) 
-9.4997(14) X 10"^'' 
33.8148(16) 
-1.35(12)xl0"^'' 
635.02377864(7) 
150.773037(2) 
8.3(14) xlO""-" 
1.47981(2) X 10"^ 
95.6398(7) 
52818.648(1) 
2210 
11.5 



64.940389248427(5) 
4.0063(6)xl0"^^ 
5.2x10^ 
5.8 
2.5 
2.6 

0.00912586(4) 
0.29 
3.11(32) 
277° 



Table 3: Proper Motions and Space Velocities of Millisecond Pulsars 



Pulsar 






Distance 


Pb 


V2D 


Rcf. 




(mas yr^ ^ ) 


(mas yr^^) 


(pc) 


(days) 


(km s"^) 




Isolated MSPs 


J0030+0451 


-5.5±0.9 


< -11 


310° 




<20 


1 


J0711-6830 


-15.55±0.08 


14.23±0.07 


860 




192±48 


2 


J1024-0719 


-34.9±0.4 


-47±1 


200 




48±13 


3 


J1730-2304 


20.27±0.06 




510 




52±13 


4 


J1744-1134 


18.804±0.015 


-9.40±0.06 


416" 




37±4 


2 


J1905+0400 


-3.80±0.18 


-7.3±0.4 


1700 




89±25 


This work 


B1937+21 


-0.46±0.02 


-0.66±0.02 


3580 




107±31 


5 


J1944+0907 


12.0±0.7 


-18±3 


1800 




197±58 


2 


J2124-3358 


-14.12±0.13 


-50.34±0.25 


322° 




87±35 


6 


J2322+2057 


-17±2 


-18±3 


790 




100±30 


7 



Binary MSPs 



,10437-4715 


121.438±0.003 


-71.438±0.007 


157° 


5.741 


146±3 


8 


J0613-0200 


1.84±0.08 


-10.6±0.2 


1700 


1.198 


100±23 


2 


J0751 + 1807 


-1.3±0.2 


-6. Oil. 8 


610 


0.263 


18±4 


9 


J1012+5307 


2.4±0.2 


-25.2±0.2 


840 


0.604 


110±25 


10 


J1023+0038 


10±1 


-16±2 


1300 


0.198 


90±26 


11 


J1045-4509 


-6.0±0.2 


5.3±0.2 


1940 


4.083 


155±36 


2 


J1455-3330 


5±6 


24±12 


530 


76.174 


57±30 


4 


J1600-3053 


-0.99±0.10 


-6.7±0.5 


2930 


14.348 


92±26 


12 


J1640+2224 


1.66±0.12 


-11.3±0.2 


1160 


175.461 


69±18 


13 


J1643-1224 


6.0±0.1 


4.1±0.4 


454° 


147.017 


22±2 


2 


J1709+2313 


-3.2±0.7 


-9.7±0.9 


1400 


22.711 


83±23 


14 


J1713+0747 


4.917±0.004 


-3.933±0.010 


1050° 


67.825 


29±3 


15 


J1853+1303 


-1.67±0.08 


-2.91±0.12 


2100 


115.654 


72±22 


This work 


B1855+09 


-2.899±0.013 


-5.45±0.02 


910 


12.327 


37±10 


16 


J1903+0327 


-2.01±0.07 


-5.20±0.12 


6400 


95.714 


180±5 


17 


J1909-3744 


-9.510±0.007 


-35.859±0.019 


1270° 


1.533 


231±13 


2 


J1910+1256 


0.21±0.10 


-7.25±0.12 


2300 


58.467 


117±34 


This work 


J1911-1114 


-6±4 


-23±13 


1220 


2.716 


129±75 


4 


J1918-0642 


-7.2±0.1 


-5.7±0.3 


1240 


10.913 


49±11 


18 


B1953+29 


~0.9±0.1 


-4.1±0.1 


4640 


117.349 


214±56 


This work 


B1957+20 


-16.0±0.5 


-25.8±0.6 


2490 


0.382 


422±110 


19 


J2019+2425 


-9.41±0.12 


-20.60±0.15 


1490 


76.512 


192±51 


20 


J2033+1734 


-5.94±0.17 


-11.0±0.3 


2000 


56.308 


162±43 


16 


J2051-0827 


5.3±1 


0.3±3 


1040 


0.099 


20±10 


21 


J2129-5721 


9.35±0.1 


-9.47±0.1 


1360 


6.625 


104±30 


2 


J2229+2643 


1±4 


-17±4 


1450 


93.016 


126±45 


22 


J2317+1439 


-1.7±1.5 


7.4±3.1 


820 


2.459 


40±15 


23 



Refere nces: flllLommen et al.1 ^20061 ) , (2)IVerbiest et al.l ^20091 ); (3)lHotan et aO l|2006l '): f 4')IToscano et al.l ||1999I ); (5)rChampion et al.l 
1.2005'); (6) 'Nice & Taylor" (1993); (7) 'Nice &: Taylor! I ll 9951) ; (8) 'Vcrbicst ot al.' (200&); (9) INicc et al." (2005'); (10) 'Langc ct al.l 
(11) Archibald et al. (2009); (12) Ord et al. (20061); (1 3) [T ohmer et al. (2005); (14) Lowandowski ct al. (2004); (15) Splavc 
1200 5); (16) Splavcr (2004); (17) .Frcirc et al.i (i201lj): (18) | Janssen et al.. (,20ia); ( 19) |Arzoumanian et al.. a994) ; (20) .Nice ct al., 
(21) (Doroshcnko et al. |200j); (22) IWolszczan et al.l | |2000| ): (23) ICamilo eTld] 11999) ! 
° Distance derived using parallax. See reference for values. 



23 



Table 4: WBMSPs used for SEP and Lorentz invariance/momentum conservation tests 



References 

p^orimer et al. (2006) 



Pulsar 


P 


Pb 


e 




h 






d 


Rcf. 




(ms) 


(days) 




n 


(Mo) 


(mas yr"*) 


(mas yr"') 


(kpc) 




J0407+1607 


25.7017 


669.0704 


9.368(6)xl0"* 


291.74(2) 


0.002893 






1.33 


1 


J0437-4715 


5.7574 


5.7410 


1.91685(5)xl0"^ 


1.22(5) 


0.001243 


121.438±0.003 


-71.438±0.007 


0.157 


2, 3, 4 


J0614-3329 


3.1487 


53.5846 


1.801(l)xl0"*' 


15.92(4) 


0.007895 






1.9 


5 


J1045-4509 


7.4742 


4.0835 


2.37(7)xl0"^ 


243(2) 


0.001765 


-6.0±0.2 


5.3±0.2 


1.94 


6 


J1125-6014 


2.6304 


8.7526 


l(13)xl0"'' 


273(87) 


0.008128 






1.50 


7 


J1216-6410 


3.5394 


4.0367 


7(59) xlO^'^ 


177(1) 


0.001669 






1.33 


7 


J1455-3330 


7.9872 


76.1745 


1.697(3)xl0"^ 


223.81(1) 


0.006272 


5±6 


24±12 


0.53 


8, 9 


J1600-3053 


3.5979 


14.3484 


1.7373(2) X 10"'' 


181.73(1) 


0.003558 


-0.99±0.10 


-6.7±0.5 


2.93 


10 


J1640+2224 


3.1633 


175.4607 


7.97262(14)xl0"^ 


' 50.7308(10) 


0.005907 


1.66±0.12 


-11.3±0.2 


1.16 


11, 12 


J1643-1224 


4.6216 


147.0174 


5.0579(4)xl0"* 


321.850(4) 


0.000783 


O.OitO.l 


4.1±0.4 


0.454 


6 


J1709+2313 


4.6312 


22.7119 


1.87(2) xlO"*^ 


24.3(6) 


0.007438 


-3.2±0.7 


-9.7±0.9 


1.4 


11, 13 


J1711-4322 


102.6183 


922.4707 


2.375(6)xl0-^ 


293.75(12) 


0.003434 






3.84 


7 


J1713+0747 


4.5701 


67.8251 


7.49406(13) xlO"" 


176.1915(10) 


0.007896 


4.917±0.004 


-3.933±0.010 


1.05 


14, 15 


J1732-5049 


5.3125 


5.2630 


9.8(20) X 10"" 


287(12) 


0.002449 






1.41 


16 


J1751-2857 


3.9149 


110.7465 


1.283(5) X 10"* 


45.52(19) 


0.003013 






1.1 


17 


J1804-2717 


9.3430 


11.1287 


3.1(5)xlO"^ 


160(4) 


0.003347 






0.78 


18, 19 


J1853+1303 


4.0918 


115.6538 


2.3686(12)xl0"^ 


346.61(4) 


0.005439 


-1.67±0.08 


-2.91±0.12 


2.1 


This work 


B1855+09 


5.3621 


12.3272 


2.170(3) X 10"^ 


276.39(4) 


0.005557 


-2.899±0.013 


-5.45±0.02 


0.91 


20, 21 


J1910+1256 


4.9836 


58.4667 


2.3017(3) X 10"'' 


106.016(9) 


0.002962 


0.21±0.10 


-7.25±0.12 


2.3 


This work 


J1918-0642 


7.6459 


10.9132 


1.991(13)xl0"^ 


218.6(4) 


0.005249 


-7.2±0.1 


-5.7±0.3 


1.24 


22 


J1933-6211 


3.5434 


12.8194 


1.3(4)xl0"'' 


116(22) 


0.011749 






0.52 


23 


B1953+29 


6.1332 


117.3491 


3.3026(5) X 10"" 


29.489(8) 


0.002417 


-0.9±0.1 


-4.1±0.1 


4.64 


This work 


J2016+1948 


64.9404 


635.0238 


1.47980(2)xl0"^ 


95.640(1) 


0.009126 


1.28±0.26 


2.82±0.34 


2.5 


This work 


J2019-5721 


3.9345 


76.5116 


1.1109(4)xl0"'' 


159.03(2) 


0.010687 


-9.41±0.12 


-20.60±0.15 


1.49 


24, 25 


,72033+1734 


5.9489 


56.308 


1.2876(6) X 10"" 


78.23(3) 


0.002776 


-5.94±0.17 


-11.0±0.3 


2.0 


21 


J2129-5721 


3.7263 


6.6255 


1.21(3)xl0"^ 


196.3(1.5) 


0.001049 


9.35±0.1 


-9.47±0.1 


1.36 


6 


J2229+2643 


2.9778 


93.0159 


2.556(2) X 10—* 


14.42(0.05) 


0.000839 


1±4 


-17±4 


1.45 


26, 27 




'l)'Lorimer ct 


al.! (120051); (2)Uohnston ( 


3t al.'('1993'); (3)'van 


Stratcn ct al.'f2001"); (4)'Vcrbiost ot al.l (12008)); (5)'Ransom et al.l(2011'): (6)'Vorbiost ct al 



(8) Toscano ot al. (1999); (9) Lorimor ct al. (1995); (10) Ord ct al. (2006); (11) 

" 7}jSta 



(2004); (14) Foster ct al. ( 1993(); (15j,Splavcr ct al., (,2005,); (16) .Edwards He Bailca (',2001.'); (1, , ^ ^^^^^^ , , 

[Scaclstcin ot al. (1986); (21) ISplaveil 120041 ); (22) IJanssen et al.l l2010f ); ('23') IJacobv et Ed] 120071 ): (24') INice et al] (|1993| )~T25') iNice 
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